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We consider elastic scattering of electrons off a proton target. The parity violating (PV) asymme- 
try arises at leading order in a due to interference of 7 and Z exchange. The radiative corrections 
to this leading mechanism were calculated in the literature and included in experimental analyses, 
except for 7Z box and cross-box contributions. We present here a dispersion calculation of these 
corrections in forward kinematics. We demonstrate that at the GeV energies of current PV experi- 
ments, such corrections are not suppressed by the small vector weak charge of the electron, as occurs 
in the atomic PV. Our results suggest that the current theoretical uncertainty in the analysis of 
the QWEAK experiment might be substantially underestimated, and more accurate account of the 
dispersion corrections are needed in order to interprete the PV data. 
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Precision tests of the Standard Model at low energies 
provide an important tool to search for New Physics and 
to constrain model parameters. Such tests involve high 
precision measurements of observables that are typically 
suppressed or precisely vanish in the Standard Model 
(SM). Prominent examples of such observables include 
the electric dipole moment and neutrino magnetic mo- 
ments. Another important example of a parameter of 
the nuclcon structure suppressed in the SM is the weak 
charge of the proton, = 1 — Asin^ 9w- With the 
value of the weak mixing angle at low momentum trans- 
fers sin^ewiO) = 0.23807 ± 0.00017 pL^, the SM predicts 
the proton weak charge of order « 0.05. A precise (4%) 
measurement of the weak charge of the proton is the aim 
of the QWEAK experiment at Jefferson Lab J^. 

In order to achieve the required precision in the 
QWEAK experiment, the radiative corrections have to 
be considered. This was done in various works, to men- 
tion the most important references [HIS], and the com- 
bined estimate of the theoretical uncertainty is currently 
2.2%. This level of precision, coupled with a 2% measure- 
ment of the parity violating asymmetry, would allow for 
a 0.3% determination of 6w at low energies. The main 
difficulty in calculating the radiative corrections origi- 
nates in the hadronic structure-dependent contributions 
from the box diagrams with the exchange of 77, ZZ, WW 
and 'yZ. Since the parity conserving amphtude at lead- 
ing order has a ^ pole, the exchange of two photons 

only leads to a correction that vanishes at = 0. This 
amounts in a contribution ~ ctQ^, with a « j^, that 
can safely be neglected. Parity violating amplitude in 
the OBE approximation has no such pole, and the re- 
spective correction remains finite in the forward direc- 
tion. The ZZ and W^W^-boxes were estimated in fl] |3] to 
give a large correction that comes from hard exchanged 
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bosons' momenta in the loop, ^ Mz{Mw), whereas low 
momenta contributions are suppressed by an extra power 
of Gi:-. In this case, all subprocesses inside the loop can 
be treated perturbatively and the contribution can be cal- 
culated reliably. The situation with jZ-hox is however 
more complex, since there is generally no reason for hard 
exchanged momenta to dominate the loop with respect 
to low momenta. For atomic PV it was observed 
that these dispersion corrections are suppressed, as the 
contributions from the box and the crossed box cancel, 
and the only non-zero term is proportional to the small 
vector charge of the electron thus leading to a correc- 
tion below 1%. In 1, 4 , this argument was adopted to 
high energy electrons, guided by the assumption of high 
momentum dominance of the loop integral. Clearly, the 
overall theory uncertainty of 2.2% relies heavily on this 
cancellation mechanism. The goal of this letter is to in- 
vestigate the dispersion correction due to 7Z-box graph 
in the kinematics of the QWEAK experiment. We will 
provide an explicit calculation of the box and crossed-box 
corrections in the framework of dispersion relations. 

Elastic scattering of massless electrons off a nucleon, 
e{k) + N{p) e{k') + N{p'), in presence of parity vio- 
lation and in most general form is described by six am- 
plitudes fi{v,t), i = 1,2,.. .,6, that depend on v = ^ 

with average momenta K = and P = and 

t = < the elastic momentum transfer, with A = 
k — k' = p' — p. The basis can be written in the following 
form: 
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where only electromagnetic and weak neutral currents are 
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considered. The amplitudes /i,2,3 are parity conserving 
(PC), and /4,5.6 are explicitly parity violating. Since we 
are interested in very forward scattering angles « 6° 
corresponding to the QWEAK kinematics [2^, the num- 
ber of structures that are relevant is further reduced. The 
magnetic terms vanish in the forward direction and can 
be neglected. The Gordon's identity allows to rewrite 
N-ft'N = j^NN - Nia^'^'^N where we made 

use of the nucleon state normalization NN — 2M. The 
parity conserving amplitude arises due to an exchange 
of at least two photons or boson, and is suppressed 
as ~ 0{t), as compared the leading PC amplitude. Fi- 
nally, amplitude /g depends on nucleon spin and makes 
no contribution to observables with an unpolarized tar- 
get. Only two amplitudes of relevance remain in the 
forward direction, and denoting their forward values as 
/,; = fi{i', t ~ 0) we obtain 

nt ^ 0) = —^2fMk')Mk) - ^Uuik')f^,u{k) (2) 

The form of the forward amplitude coincides with that 
for spin-0 target [5j. To leading order in t and in Fermi 
constant, the parity violating asymmetry arises from the 
interference of PC and PV amplitudes, 
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At tree level, the amplitudes reduce to /f = 1 and 
fOBE ^ ge^QP^ ^ I - 4sin2 Ow, thus the PV asymmetry 
gives a direct access to the weak charge of the proton. 
Using the exact forward values of the amplitudes (the 
neglect of the t-dependence of respective form factors) 
only leads to power corrections in t and can be safely 
neglected [T]. In the following, we will assume that all 
the radiative corrections except the 7Z direct and crossed 
boxes are known. Denoting these other corrections by 
SrCi we introduce the correction due to 7Z exchange as 



37Z 



f fOBE 
fOBE 



^RC, 



with d-yz a complex function energy v. Its real part con- 
tributes to the parity violation asymmetry as 



47ra\/2 



QPy[l + ReSRc + ReS^z{iy)] (5) 



To calculate the real part of the 7Z direct and crossed 
box graphs shown in Fig. [T| we adopt dispersion rela- 





FIG. 1: Direct and crossed box diagrams 
tion formalism, and we start with the calculation of the 



imaginary part of the direct box (the crossed box contri- 
bution to the real part will be calculated using crossing 
behavior) , 



ImT, 



(6) 



where Q'^ — — (fc — fci)-^ denotes the virtuality of the 
exchanged photon and Z (in the forward direction they 
carry exactly the same Q^), and we explicitly set the in- 
termediate electron on-shell. In the cm. of the (initial) 
electron and proton, Ei = 1 with s the full cm. 

energy squared and W the invariant mass of the interme- 
diate hadronic state. Note that for on-shell intermediate 
states, the exchanged bosons are always spacelike. The 
leptonic tensor is given by 



In 



(7) 



In the case of the elastic hadronic intermediate state, 
the needed structure A{e) x V{p) always contains the 
explicit factor of or gy. Correspondingly, the cor- 
rection 5~^z is not enhanced with respect to the small 
tree-level coupling, and is generally small, in accordance 
with [Ij. We therefore turn to the inelastic contribution. 
In the forward direction, the imaginary part of the doubly 
virtual "Compton scattering" (7*^ Z*p) amplitude is 
given in terms of the structure functions ^1,2,3(2;, Q^), 
with X = the Bjorken variable. Making use of gauge 
invariance of the leptonic tensor, we have 

W^- ^ Jd^xe^'^^ <p\T{j:l{x)rz{0)}\p> (8) 

pp. pu _ 
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Contracting the two tensors, one obtains after little 
algebra 
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with = (M m^)2 the pion production threshold, 
and the Q^-integration is constrained below a maximum 
value Qmax = (^-M )(s-w ) ^ ^ condition of on-shell 
intermediate states for an imaginary part calculation. In 
order to write down dispersion relation for the function 
S'rz{'^)i we consider its behavior under crossing. Crossing 
corresponds to a CP-transformation applied to a part 
of the amplitude, so that it relates the original reac- 
tion e^(fc) + N{p) — > e~(fc') + N{p') to the reaction 



e+i-k') + N{p) 



(— fc) + N(p'). The requirement 



that the crossed reaction be described by the same invari- 
ant amplitudes taken at the crossed kinematics v ^ —v 
{K —t ~K with P unchanged), imposes constraint on 
the form of its zz-dependence. The tensor that multiplies 
the amplitude is even under crossing (being an axial 
vector). Under C-parity applied to the electron part, the 
tree level contribution to is also even, and as func- 
tion of V the OBE amplitude is an even function - as it is 
observed, in fact, since it only depends on the elastic mo- 
mentum transfer that is unchanged under crossing. At 
one-loop order, exchange of two vector currents (e.-m. 
and vector part of NC) leads to C-even behavior, and 
vector-axial vector exchange to C-odd. Correspondingly, 
the part of ImS^z that contains g\ is an odd function of 
v, whereas the one with gy is even. We distinguish then 
the two contributions 5^Zv 5^Za that obey dispersion 
relations of two different forms. 
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This is just another formulation of the mechanism of 
the cancellation between the box and crossed box graphs 
observed in [S]. For small values of that are rele- 
vant for PV in atoms, the explicit factor of v in front 
suppresses RbS^Za- Instead, no cancellation occurs for 
Red^Zy which is however suppressed by the small value 



of = -1 4sin2 ew The result of Eq. (flO I can be 



considered as a sum rule, since it represents the quan- 
tity that is to be measured (the proton's weak charge 
plus corrections to it) through other observables (PVDIS 
structure functions), and this relation does not rely on 
any assumption, other than the neglect of higher order 
radiative corrections. 

In absence of any detailed experimental data on PVDIS 
structure functions, to provide estimates of Re(5^2^(zy), 
we need a phenomenological model for and de- 
pendence of Fi^2- Since we need input from low to 
high values in both variables, one has to distinguish 
two different contributions: nucleon resonances contri- 
butions in the intermediate states, and high energy non- 
resonant part. It was shown in a number of papers that 
a combined color dipole picture/ vector meson dominance 
(CDP/VMD) approach allows to describe world data on 
the parity-conserving DIS structure functions Fi^2 over a 
wide range. Although no data is available to directly 
judge, whether or not the PVDIS structure functions 
Fi^2 would follow exactly the same pattern, we can ex- 
pect that they may be similar. In CDP, the interac- 
tion of the high-energetic photon with the proton target 
occurs in two steps. Firstly, the photon fluctuates into 
a quark-antiquark pair that forms a color dipole. Sec- 
ondly, this dipole interacts with the target by means of 
exchange of gluons, as shown in Fig. [2j Correspondingly, 
the color dipole-proton scattering cross section is the 
only universal non-perturbative ingredient, and all the 




FIG. 2: Color dipole picture 



information about the external particles is contained in 
the hadronic (qq) wave functions of the virtual photon 
or vector meson, depending on the reaction under con- 
sideration. In our case, it is the wave function of the 
Z-boson, and its vector part is the same as for a vir- 
tual photon of equal virtuality (which is the case due 
to forward kinematics), except for the quark flavor de- 
pendence. This flavor dependence goes as for the 
77 DIS structure functions, and as J2q ^q9v ^'^^^ 
PVDIS, provided we restrict ourselves to the lightest fla- 
vors (since all intermediate states have to be real parti- 
cles, production of a heavy quark-antiquark pair requires 
very high energy) and the dependence on quark masses 
can be safely neglected in the photon and Z wave func- 
tions. By taking the SM values, J2q=u ds^q^ whereas 

E,=n,.;,.egff?.= i(l + (l-4sin2 0H'))« I We also re- 
mind the reader that the very first PVDIS measurment 
on the deuteron target [3] was used to confirm the par- 
ton model: in the scaling limit, the isoscalar structure 
functions of DIS and PVDIS depend on the same combi- 
nation of the quark PDFs that cancel out in the asym- 
metry, leaving only kinematical factors and couplings. 

This similarity is the case, as well, for the A (1232) res- 
onance, as shown in Ref. [7]. Following a simple isospin 
decomposition, one obtains 



<^\Jnc\p> 



< A|(2-4sin2 ew)JL 



(11) 



with / the strong isospin (u,d-quarks are only kept). 
Since A has isospin 3/2 and the proton 1/2, the isoscalar 
correction is small, and the NC excitation of the A 
only differs from the electromagnetic one by a factor of 
~ 1.08. 

We conclude that the assumption of the similarity of 
the interference 7Z structure functions to the usual ones 
with two virtual photons is supported within the parton 
model and CDP at high energy, and at low energies, at 
least for the most important A(1232) resonance. We will 
use this assumption to provide a realistic estimate for the 
dispersion correction to the proton weak charge. Alter- 
natively to the DIS structure functions, one can use the 
transverse and longitudinal virtual photon cross sections 
aT,L related to Fi^2 as (see, e.g., Ref. [9^ for details) 
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F, 
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To proceed, we match the relative strength of the res- 
onances in Breit-Wigner form and the high energy Regge 
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part to fit the real Compton data (for the parameters see 

a,,{W ,0) = l.(^2„M2)2 + M2r|+^-^ ) 

(13) 

Further, we employ two different forms of depen- 
dence for the two contributions. Ref. |10j gives an- 
alytical form of W"^ and dependence of the longi- 
tudinal and transverse virtual photon cross sections as 

ing variable ry = l^i'{w^ ) ^^"^ '^o = "niQ^ — 0); so that 
for the real photons the known Regge asymptotics is re- 
covered, and we refer the reader to Ref. [TU] for fur- 
ther details. For the resonances, transition form fac- 
tors are used. The latter are to some extent known for 
a number of resonances and we assume a dipole form, 

FTiQ"^) = (i+q1/a^)^ ' and i^L(Q^) = (i+QVA^)^.5 , with 
Awl GeV. Finally, the form that is used for numerical 
estimates is 

-I- a^Js'^(M^2 (^4) 

We present results of the dispersion calculation in Fig. 
|3] It can be seen that starting from Eiab ~ 1 GeV, 
high energy (Regge) contribution dominates the contri- 
bution from the resonances. This is the consequence of 
a relatively slow convergence of the dispersion integral 
for Regge part, while the resonances drop very fast. In 
the presented calculation, the upper limit of the inte- 
gration over v' was chosen to be 500 GeV, although the 
weighting ensures the convergence already at lower 
values. While at very low energies the correction is in- 
deed very small, at the 1.16 GeV energy of the QWEAK 
experiment the correction is 5.7%. More specifically, 
QWEAK aims at comparing the measured weak charge 
of the proton, iro^^py ^jg^^ie as given in SM, 

[1 -t- (5flC + ReJ-yz] and from this comparison draw 
conclusions about the New Physics contributions. The 



current estimate of the uncertainty due the the correc- 
tions in the square brackets is 2.2%, and this estimate 
relies on the assumption that S^^z is highly suppressed 
(<0.65%). As explained above, this estimate is taken 
over from low energy estimates for PV in atoms, and is 
not based on any microscopic calculation. Although the 
numbers presented here are themselves model-dependent, 
our calculation shows that the 7Z box diagrams can be 
almost an order of magnitude larger than it was believed 
to date, and this result suggests larger possible theoreti- 
cal errors for the QWEAK experiment. If the uncertainty 
in the dispersion correction is to be comparable to the 
proposed 2% experimental error in one may need 

to calculate the dispersion 7Z correction (that we think is 
near 6%) to a fractional accuracy of order 30%. Alterna- 
tively, uncertainties in these dispersion corrections could 
provide a limit on the precision of a Standard Model test. 
Since the calculation uses the PVDIS structure functions 
as input, it would be extremely helpful to have experi- 
mental data on PVDIS to check the model adopted here. 




(GeV) 

FIG. 3: Results for ReJ^z^ as function of energy. The con- 
tributions of nucleon resonances (dashed line), Regge (dash- 
dotted line) and the sum of the two (solid line) are shown. 
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